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Abstract. We report on determinations of the running mass for charm quarks from deep-inelastic 
C^ • scattering reactions. The method provides complementary information on this fundamental param- 

CsJ ' eter from hadronic processes with space-like kinematics. The obtained values are consistent with 

but systematically lower than the world average as published by the PDG. We also address the 

consequences of the running mass scheme for heavy-quark parton distributions in global fits to 

deep-inelastic scattering data. 
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Quark masses are fundamental parameters of the gauge theory of the strong interac- 
5J^| tions, Quantum Chromodynamics (QCD). They are, however, not directly observable 

7^, ' due to confinement. Rather one has to compute the dependence of cross sections or 

other measurable quantities on the heavy quark mass including higher order radiative 
corrections and renormalization - a procedure which requires the choice of a scheme 

Q^ ■ for the definition the mass parameter. There exists a variety of schemes for heavy-quark 

^ ! masses. The most popular ones are the on-shell and the MS scheme. In the former the 

so-called pole-mass m coincides with the pole of the heavy-quark propagator at each 
order in perturbative QCD. This definition, however, has intrinsic theoretical limitations 

Q . with ambiguities of order (^(Aqcd) and a strong dependence of the value of the mass 

parameter on the order of perturbation theory. The MS scheme is one of the so-called 
short-distance masses for heavy quarks. It realizes the concept of a running mass m{iJ.) 
depending on the scale jl of the hard scattering in complete analogy to the treatment of 

^^ . the running coupling aj(ju). As a benefit, predictions for hard scattering cross sections 

H I in terms of the MS mass display better convergence properties and greater perturbative 

stability at higher orders. 

Cross sections for the production of heavy-quarks in deep-inelastic scattering (DIS) 
are particularly well suited to confront the quark mass dependence of theoretical pre- 
dictions in perturbative QCD with experimental measurements in space-like kinematics. 
For the production of charm quarks in neutral (NC) or charged current (CC) DIS there 
exists precise data from the HERA collider and fixed-target experiments. In QCD the 
DIS heavy-quark structure functions Fj^ which parametrize the hadronic cross section 
are subject to the standard factorization 



Fk{x,Q\m^)= "£ 



/KM^)®Q.■(e^«^av(M^)) 



(x), ;t= 1,2,3, (1) 
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where the perturbative coefficient functions Q,/ are known to next-to-leading order 



(NLO) for CC [1,2] and approximately to next-to-next- to-leading order (NNLO) for 
NC [3, 4]. Q^ and x are the usual DIS kinematical variables and m is the heavy-quark 
(pole) mass. In eq. (1) we also display all dependence on the other non-perturbative 
parameters, i.e. the parton distribution functions /, (PDFs) for light quarks and gluons as 
well as the strong coupling constant Us. The running mass definition can be implemented 
in eq. (1) to the respective order in perturbation theory by simply following the standard 
procedure for changing the renormalization condition, i.e. m — )■ m{\x) see ref. [5]. 

The present study has two aspects: First of all we determine the MS charm mass 
mc(mc) from DIS in a variant of the global analysis [6] and compare to the world average 
as published by the PDG. Secondly, we investigate the improvements in the uncertainty 
of heavy-quark PDFs in a global fit within the fixed-flavor number scheme (FFNS) when 
the running mass scheme is applied. 

1. Running charm-quark mass 

The parametric dependence of the DIS structure functions F^ in eq. (1) on m can 
be used for a determination of the heavy-quark mass. The sensitivity of this proce- 
dure relates directly to the corresponding uncertainty on the measurements of F^. E.g., 
for charm production in NC DIS the nucleon structure function F2 yields Amc/nic ~ 
O.75AF2/F2, which implies that an experimental accuracy of 8% for F2 translates into 
an uncertainty of 6% for the charm-quark mass [5]. With the precision of current DIS 
data for charm production this suggests an error on mc{mc) of (^(few)% as the ultimate 
precision in the approach based on inclusive structure functions. 

Starting from eq. (1) we have extracted the MS charm mass mc{mc) in a phenomeno- 
logical study similar to [7], i.e. a global fit including fixed-target (CCFR [8], NuTeV [9]) 
and collider data [10, 1 1] in the FFNS (with Uf = 3) as a variant of the ABKM one [6]. 
We have obtained [5], 

mcinic) = 1.26 ±0.09(exp) ± 0.11 (th)GeV at NLO, (2) 

mcinic) = 1.01 ±0.09 (exp) ±0.03 (th)GeV at NNLOapprox , (3) 

to NLO and approximate NNLO in perturbation theory. The quoted experimental un- 
certainty results from the propagation of the statistical and systematic errors in the data 
with account of error correlations whenever available. The theoretical uncertainty is es- 
timated from the variation of the renormalization and factorization scales for the choice 
/i^^ = jlr^ = Q^ + Kin^. for F2 with K in the range K G [0, 8] . For consistency it has been 
checked that different scale choices do not deteriorate the statistical quality of the fit. Be- 
sides the charm mass nic quoted in eqs. (2) and (3) and the value of the strong coupling 
constant determined to as(Mz) = 0.1 135 ±0.0014 at NNLO [6], the global fit relies on 
the same 22 PDF parameters as in ABKM [6]. This is particularly important in order 
to control simultaneously the sensitivity to the strange-quark PDF in CC case which is 
correlated with ntc due to the Bom process W^s -^ c. 

At this conference, NOMAD has reported a new analysis of CC fixed-target data. The 
increased precision due to the very high statistics has been used to determine the running 
mass mc{mc) to NLO exclusively from CC DIS leading to the (preliminary) result [12] 

mc{mc) = 1.070 ±0.067 (exp) ±0.050(th)GeV at NLO . (4) 



One has to compare the numbers in eqs. (2), (3) and (4) with the world average of 
the PDG [13] quoted in the MS scheme as mcinic) = 1.27 +qq^ GeV, which is entirely 
based on lattice computations or analyses of experimental data with time-like kinematics 
from e+e -collisions, e.g. with the help of QCD sum rules. It is therefore interesting to 
note that the DIS results in eqs. (2), (3) and (4) for hadronic processes with space-like 
kinematics are consistent with but systematically lower than the PDG value. In order 
to understand at least one source of this deviation, one should note the QCD sum rules 
analyses typically assume the Bethke world average for the value of the strong coupling 
constant [14], which is as{Mz) = 0. 11 84 ± 0.0007 and therefore significantly larger 
than the ABKM result, see also [15] and references therein. A recent determination 
from a charmonium QCD sum rules analysis [16] quotes a central value mc{mc) = 
1.257 ±0.026 GeV using the world average and parametrizes separately the dependence 
of mc{mc) on value of as{Mz). Using the ABKM value as{Mz) = 0.1135 instead one 
extracts from Table 15 in [16] a value of mc{mc) = 1.257 ±0.026 GeV, i.e. a systematic 
shift downwards at the level on Ic statistical uncertainty bringing QCD sum rule 
analysis from e+e^-coUisions in better agreement with the DIS results in eqs. (2), (3) and 
(4). Note that the latter determinations account for the full correlation of the dependence 
on nic and a^ through a simultaneous fit of the parameters. 

Including the PDG constraint on mc{mc) into the fit of ref. [5] we get the value of 
mc{mc) = 1.18 ± 0.06 GeV at NNLO. The recent data on the charm structure function 
^2*^ at low Q^ published by the HI collaboration [17], which were not included into 
the fit of ref. [5], are compared with the running-mass scheme predictions of ref. [5] 
obtained with this value of mc{mc) in fig. 1 (left). The agreement is quite good therefore 
these data also prefer somewhat smaller value of mc{mc) than the PDG average. 

2. Heavy quark PDFs 

At asymptotically large scales Q ^ mc,mt the genuine heavy quark contributions 
in a FFNS with nf = 3 grow as as{Q^)ln{Q^ /m^) and can be resummed by means 
of standard renormalization group methods. This procedure leads to so-called heavy 
quark PDFs in theories with effectively nf = 4 and nf = 5 light flavors, which are 
the appropriate descriptions for processes at LHC energies. The PDFs for charm- and 
bottom-quarks in nf = 4- and nf = 5-flavor schemes are generated from the light flavor 
PDFs in a nf = 3-flavor FFNS as convolutions with massive operator matrix elements 
(OMEs), see e.g. [6]. The uncertainty on heavy-quark PDFs is therefore directly related 
to the accuracy of the quark masses m^ or m^, which appear parametrically in the OMEs. 
This uncertainty can be significantly reduced through the use of the MS scheme, which, 
of course, has to be applied also to the massive OMEs. 

In our global fit [5] mc{mc) has been left a free parameter supplemented by the PDG 
constraint. In this way we have generated a charm-PDF with comparable uncertainties 
to the one of [6] (which has used the pole mass definition for nic). For m^ the currently 
available DIS data displays no sensitivity at all and we have constrained mi,{mi,) directly 
to its PDG value [13], i.e. mb{mb) = 4.19+00^ GeV. As shown in fig. 1 (right) the 
uncertainty of the resulting bottom-PDF is greatly reduced. This improvement will 
certainly have an impact on LHC phenomenology, e.g. allowing for precise predictions 
for the production of single-top-quarks. 
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FIGURE 1. Left: The data on the charm structure function F2" by the HI collaboration [17] in compar- 
ison with the running-mass scheme predictions of ref. [5]. Right: The /j-quark PDF uncertainties obtained 
in the global fits: The dotted (red) lines denote the ±la band of relative uncertainties (in percent) and 
the solid (red) line indicates the difference in the central prediction resulting from the change of the mass 
scheme and using mh{mh) — 4.19 GeV (bottom). For comparison the shaded (grey) area represents the 
uncertainties in the ABKM fit [6]. 



REFERENCES 



1. 

2. 
3. 
4. 
5. 

6. 



9. 
10. 
11. 
12. 

13. 
14. 
15. 
16. 

17. 



T Gottschalk, Phys.Rev. D23, 56 (1981). 

M. Gluck, S. Kretzer, and E. Reya, Phys.Lett. B380, 171-176 (1996), hep-ph/ 9 603304. 

E. Laenen, S. Riemersma, J. Smith, and W. L. van Neerven, Nucl. Phys. B392, 162-228 (1993). 

N. LoPresti, H. Kawamura, S. Moch, and A. Vogt(2010), arXiv: 1008 .0 951. 

S. Alekhin, and S. Moch, Phys.Lett. B699, 345-353 (201 1), arXiv :1011.5790. 

S. Alekhin, J. Blumlein, S. Klein, and S. Moch, Phys.Rev. D81, 014032 (2010), 

arXiv:0908.2766. 

S. Alekhin, and S. Moch, Phys. Lett. B672, 166-171 (2009), arXiv: 0811 . 1412. 

A. Bazarko, et al., Z.P/!yi. C65, 189-198 (1995), hep-ex/9406007. 

M. Goncharov, et al., P/iys.TJev. D64, 112006 (2001), hep-ex/0102049. 

S. Chekanov, et al., Phys.Rev. D69, 012004 (2004), hep-ex/0308068. 

F D. Aaron, et al., Eur. Phys. J. C65, 89-109 (2010), arXiv: 907 .2 64 3. 

R. Petti (201 1), these proceedings, 

https : //wiki . bnl . gov/ conferences /index . php/DIS-2 011 

K. Nakamura, et al., J.Phys.G G37, 075021 (2010). 

S. Bethke, Eur.Phys.J. C64, 689-703 (2009), arXiv: 0908 . 1135. 

S. Alekhin, J. Blumlein, H. Bottcher, and S. Moch (201 1), arXiv :1104.0469. 

B. Dehnadi, A. H. Hoang, V. Mateu, and S. Zebarjad (2011), arXiv: 1102 .22 64. 
F Aaron, et al. (201 1 ), a r Xi v : 1 1 6 . 1 2 8 . 



